Aging in coherent noise models and natural time 
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Event correlation between aftershocks in the coherent noise model is studied by making 
use of natural time, which has recently been introduced in complex time-series analysis. 
It is found that the aging phenomenon and the associated scaling property discovered 
in the observed seismic data are well reproduced by the model. It is also found that 
the scaling function is given by the g-exponential function appearing in nonextensive 
statistical mechanics, showing power-law decay of event correlation in natural time. 
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In recent years, there has been an increasing interest in extended dynamical systems exhibiting 
avalanches of activity, whose size distribution is scale-free. Examples of such systems are earthquakes [1], 
, rice piles [2], extinction in biology [3], evolving complex networks [4], and so on. Up to now, there is no 
c/3 unique and unified theory for such systems, but one of the candidates may be the notion of self-organize 
criticality (SOC) [5]. A key feature common in SOC models is that the whole system is under the influ- 
J2 ' ence of a small driving force that acts locally. These systems evolve towards a critical stationary state 
having no characteristic spatio-temporal scales without fine-tuning parameters. Two celebrated examples 
' of such extended driven dynamical systems are the Olami-Feder-Christensen model of earthquakes [6] and 
y _ the Bak-Sneppen model for biological evolution [7]. 

Q . On the other hand, there exists another kind of simple and robust mechanism producing scale-free 
^ ' behavior in the absence of criticality. An important example is the coherent noise model [8], which 
we shall study here. The coherent noise model has been introduced to describe large-scale events in 
^ ' evolution. It is based on the notion of external stress coherently imposed on all agents of the system 
^ , under consideration. Since this model does not contain any direct interaction among agents, it does not 
QQ ' exhibit dynamical criticality. Nevertheless, it yields a power- law distribution of event (i.e., avalanche) 
0^ ', size s, which is defined by the number of agents that change their states at each time step. The model 
. allows existence of aftershocks. This is a direct consequence of the fact that, in the coherent noise model, 
1^ ' the probability of a large event to occur is increased immediately after a previous large event. 

The coherent noise model is defined in a rather simple manner. Consider a system, which consists of 
N agents. Each agent i has a threshold Xi against external stress rj. The threshold levels are chosen at 
random due to some probability distribution pthresh[x). The external stress is also chosen randomly due 
' to another distribution psjress(?7)- An agent becomes eliminated if it is subjected to the stress rj exceeding 
\ the threshold for the agent. In practice, dynamics of the model can be summarized as the following three 
I ■ steps: (i) at each time step, a random stress 77 is generated from pstress iv) ^^'^ the agents with Xi < rj 
''O ' become eliminated and are replaced by new agents with new thresholds drawn from pthresh{x), (ii) a 
Ch , small fraction / of the N agents should be chosen at random and given new thresholds, and then (iii) go 
back to (i) for the next time step. Here, (ii) corresponds to the probability for the / fraction of the whole 
agents of undergoing spontaneous transition. This is necessary for preventing the model from grinding 
to a halt (see e.g. Ref. [9] for details). 

The coherent noise model has been applied to the problems of earthquakes [9,10]. There, the Gutenberg- 
Richter law for the relation between frequency of the events and the values of the magnitude [11] and 
' the Omori law for temporal decay pattern of aftershocks [12] are found to be realized. In particular, 
the Omori law, which is of relevance to our subsequent discussion, states that the rate of the number of 
aftershocks, p, after a mainshock at t = obeys 
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Empirically, the exponent r takes a value between 0.6 and 1.5. This describes the slow power-law decay, 
and each of the relevant parts of the seismic time series is nonstationary. Such a time interval is referred 
to as the Omori regime. 

Quite recently, the physical properties of correlation in the seismic time scries have been studied in [13] 
based on analysis of the observed seismic data. It was shown that the aging phenomenon occurs inside 
the Omori regime but it disappears outside. The definite scaling property has also been identified for 
the event correlation function. These results suggest that nonstationarity of the Omori regime may have 
something to do with glassy dynamics. 

A point of crucial importance regarding the aging phenomenon of earthquake aftershocks is that, unlike 
ordinary discussions of the aging phenomenon, the two-point correlation function is defined in the domain 
of "natural time" [14], not conventional time. The concept of natural time is a kind of an internal clock 
coimting the discrete event number. It has successfully been applied to revealing physical essence of 
complex time series such as seismic electric signals, ionic current fluctuations in membrane channels and 
so on [14]. The fact that natural time is more fundamental than conventional continuous time is still 
empirical. In this respect, one may recall that the concept of continuum is recognized through continuous 
physical processes, whereas in the case of earthquakes one is concerned with a series of discrete events. 
However, clearly, more investigations are needed for deeper understanding of natural time. 

Here, wc discuss the physical properties of the Omori regime in the coherent noise model. We shall 
see that aging and scaling discovered in the observed seismic data are well reproduced by the model. 
The correlation function is found to have the form of the g-exponential function (see below), showing 
slow power-law decay. These results have also striking similarities with those recently obtained for 
the nonextensive Hamiltonian system of an infinite-range coupled rotors in the course of nonequilibrium 
relaxation dynamics [15] in conventional time, although the systems as well as the chosen random variables 
are completely different from each other. This, in turn, indicates that the aging/scaling phenomena may 
be quite universal in extended dynamical systems out of equilibria, no matter if critical or not. 

We have carried out the numerical study of the aging phenomenon in the coherent noise model with 
the exponential distribution for the external stress 

Pstressin) = exp (^—^ (a > 0) , (2) 

and the uniform distribution pthresh{x) (0 < a; < 1) for the threshold level. The results obtained are in 
order. 

First of all, we present in Fig. 1 a typical subinterval of the obtained time series of activity, where 
aftershocks are clearly identified. In Pig. 2, the rate of the probability of finding aftershocks larger than 
si following the initial large event at f = is plotted with respect to the elapsed time t. The straight 
line represents the Omori law, which allows us to identify the Omori regime. We have ascertained that 
this result is insensitive to the threshold value s\. 

To investigate the property of event correlation, following the idea proposed in [13], we have employed 
as the basic random variable the time of the nth aftershock with an arbitrary avalanche size, i„, where n 
is aforementioned natural time in the setting of our problem. The two-point correlation function is given 

by 

C{n + n^,n^) = ^*"+"-^"-^ " , (3) 

where the average is understood as the ensemble average taken in association with a number of numerical 
runs, and the variances in the denominator are given by, = (t^) — {tm}^ ■ It is noticed that this 
averaging procedure is different from that in [13], in which the time average is used. Since the Omori 
regime is nonstationary, the correlation fimction depends not only on n but also on In Fig. 3, we 
present the plots of C {n + nw,n.uj) for several values of "natural waiting time" n^. There, the clear aging 
phenomenon can be appreciated. Furthermore, as shown in Fig. 4, collapse of these curves can nicely be 
realized, following the scaling relation 

C{n + n^,n^) = C (^] , (4) 
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where C is a scaling function and a is numerically a ~ 1.05. This form of scaling is remarkably similar 
to the one presented in [13], where the observed seismic data are employed. 

To find the form of the scaling function, we have examined the semi-g-log plot of the curve in Fig. 4, 
where the g-logarithmic function is defined by 

ln,(x) = ^^^^ > 0) , (5) 

which is the inverse function of the g-exponential function 

e,{x) = [l + {l-q)xf_l^'-'^^ (6) 

with the notation [z\+ = max{0, 2}. (In the limit 5 — > 1, Inq(x) and eq{x) converge to the ordinary 
logarithmic and exponential functions, respectively.) This pair of functions is known to play a central 
role in nonextensive statistical mechanics [16,17]. The straight line of the correlation function in the scmi- 
g'-log plot shown in Fig. 5 indicates that the scaling function is given by the gr-exponential. Therefore, 
event correlation decays slowly, following a power law. This reminds us of the recent discussion [15] 
about the nonextensive Hamiltonian system in the course of nonequilibrium relaxation, though the time 
variable employed there is the conventional one. 

In conclusion, we have studied the physical properties of event correlation in the Omori regime of 
aftershocks in the coherent noise model. We have found that aging and scaling in natural time discovered 
in the observed seismic data can be reproduced remarkably well by this model. We have also found that 
the scaling function is given by the g-exponential function, and thus event correlation decays according 
to a power law. 
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Figure Captions 



Figure 1: A section of tlie time series of activity in the dimensionless units. 5000 time steps are shown 
after the largest shock out of the time series of 30000000 events. The aftershocks following the mainshock 
at t = are clearly recognized. 

Figure 2: A histogram of the time distribution (i.e., the rate) of 90000 aftershocks larger than s\ = 10 

following the mainshock at i = 0. The strength of the external stress in cq.(2) is a = 0.001, and the 
fraction / = 5 x 10~^ of iV = 200000 agents is chosen at random. It obeys a power law in eq.(l) with 
T ~ 1. All quantities are dimensionless. 

Figure 3: Dependence of the event correlation function C (n + n^i, /J«i) of the aftershocks larger than 
Si = 1 on natural time. The strength of the external stress in eq.(2) is a = 0.2 and the fraction / = 0.01 
of iV = 10000 agents is chosen at random. The ensemble average over 120000 realizations is performed. 
The values of natural waiting time are = 250, 500, 1000, 2000 and 5000 from bottom to top. All 
quantities are dimensionless. 

Figure 4: Data collapse for the correlation function C {n + nw,nui) shown in Fig. 3. The gray solid 
line corresponds to [—0.7n/nl;^^) with q ~ 2.98. 

Figure 5: The semi-g-log plot of the collapsed curve with = 1000, 2000 and 5000 in Fig. 4. The 
straight line shows that the scaling function is of the g-exponential form with q ~ 2.98. 
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